Significant Conditions on the Two-electron Reduced Density Matrix from the 
Constructive Solution of iV-representability 
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We recently presented a constructive solution to the iV-representability problem of the two- 
electron reduced density matrix (2- RDM) — a systematic approach to constructing complete condi- 
tions to ensure that the 2-RDM represents a realistic iV-electron quantum system [D. A. Mazziotti, 
Phys. Rev. Lett. 108, 263002 (2012)]. In this paper we provide additional details and derive 
further JV-representability conditions on the 2-RDM that follow from the constructive solution. The 
resulting conditions can be classified into a hierarchy of constraints, known as the (2, g)-positivity 
conditions where the q indicates their derivation from the nonnegativity of g-body operators. In 
addition to the known Tl and T2 conditions, we derive a new class of (2,3)-positivity conditions. 
We also derive 3 classes of (2,4)-positivity conditions, 6 classes of (2,5)-positivity conditions, and 
24 classes of (2,6)-positivity conditions. The constraints obtained can be divided into two general 
types: (i) lifting conditions, that is conditions which arise from lifting lower (2, g)-positivity con- 
ditions to higher (2,5+ l)-positivity conditions and (ii) pure conditions, that is conditions which 
cannot be derived from a simple lifting of the lower conditions. All of the lifting conditions and the 
pure (2, £jr)-positivity conditions for q > 3 require tensor decompositions of the coefficients in the 
model Hamiltonians. Subsets of the new TV-representability conditions can be employed with the 
previously known conditions to achieve polynomially scaling calculations of ground-state energies 
and 2-RDMs of many-electron quantum systems even in the presence of strong electron correlation. 

PACS numbers: 31.10.+Z 



I. INTRODUCTION 

Because electrons are indistinguishable with pairwise 
Coulomb interactions, the energies and properties of 
many-electron atoms and molecules can be evaluated 
from a knowledge of the two-electron reduced density 
matrix (2-RDM) [H-Q- Minimizing the ground-state en- 
ergy as a functional of the 2-RDM, however, requires 
non-trivial constraints on the 2-RDM to ensure that it 
represents an TV-electron system (N- represent ability con- 
ditions) [H-tzjj. While advances in theory and computa- 
tion enabled the accurate variational calculation of the 2- 
RDM for a variety of strongly correlated systems in chem- 
istry and physics from polyaromatic hydrocarbons [25], 
I2H ] to quantum dots [27]], the known iV-representability 
conditions for the 2-RDM, albeit rigorous, remained in- 
complete. Recently, we presented a constructive solu- 
tion to the iV-representability problem — a systematic ap- 
proach to constructing complete A^-representability con- 
ditions on the two-electron reduced density matrix (2- 
RDM) — as well as examples of new A^-representability 
conditions 28]. In the present paper we present addi- 
tional details as well as further conditions on the 2-RDM 
that follow from the constructive solution. 

The advantage of reduced variables such as the 2-RDM 
and the one-electron density is that, unlike the wavefunc- 
tion expanded in terms of determinants, their degrees of 
freedom grow polynomially with the size of the quantum 
system [3| even when the electrons are strongly corre- 
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lated [29|,|30|. Direct calculation of the reduced variables, 
however, requires that they and their functionals be con- 
sistent with a realistic A^-electron quantum system; in 
other words, the reduced variables and functionals must 
be representable by the integration of an A^-electron den- 
sity matrix. Such consistency relations are known as the 
N-representability conditions [i"14l8l. [2fj|-|24|. These con- 
ditions are particularly important to 2-RDM methods 
where they enable the direct calculation of the 2-RDM 
without the wavefunction, but they are also implicit in 
the design of realistic approximations to the density func- 
tional in density functional theory [3l|, HH . 

Minimizing the many-electron energy as a functional 
of the 2-RDM without 7V-representability conditions pro- 
duces an energy that is much lower than the exact 
ground-state energy of the quantum system. The energy 
is too low because both the energy and the computed 
2-RDM are not realistic — they are not iV-representablc. 
In the early 1960s the search for the set of necessary and 
sufficient iV-representability conditions became known as 
the iV-representability problem [H . Three important con- 
straints, known as the D, Q, and G (or 2-positivity) con- 
ditions, were developed by Coleman [4| and Garrod and 
Percus Q. The D, Q, and G conditions restrict the prob- 
ability distributions of two electrons, two holes (where 
a hole is the absence of an electron), and an electron- 
hole pair to be nonnegative. Each condition can be ex- 
pressed in the form of constraining a matrix to be positive 
semidcfinite. A matrix is positive semidefinite if and only 
if its eigenvalues are nonnegative. 

In 1978 Erdahl [8( discovered two additional semidef- 
inite constraints on the 2-RDM known as the Tl and 
T2 (or partial 3-positivity) conditions [H, [U El SI], 
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which are derivable from the nonnegativity of the three- 
electron probability distributions. Finally, Weinhold and 
Wilson [3iJ , Yoseloff and Kuhn [35[ , McRae and David- 
son (36j . and Erdahl Q derived necessary conditions on 
the diagonal part of the 2-RDM. These diagonal condi- 
tions were shown, in the context of the Boole optimiza- 
tion problem [371 ] . to be part of a complete set of clas- 
sical N -representability conditions on the two-electron 
reduced density function which is the dia gon al part of 
the 2-RDM in a coordinate representation [38j. Despite 
the solution of the classical problem, the complete set 
of quantum 7V-representability conditions remained un- 
known except for the D, Q, G, Tl, and generalized T2 
conditions as well as unitary transformations of the clas- 
sical iV-representability conditions. In 2001 Mazziotti 
and Erdahl [ill ] presented a systematic generalization of 
these constraints known as the p-positivity conditions and 
in 2002 Mazziotti [13, 39] introduced the lifting condi- 
tions; however, except for the conditions given above, 
the p-positivity conditions and the lifting conditions de- 
pend upon not only the 2-RDM but also higher-particle 
RDMs. 

The constructive solution to the A-representability 
problem provides a systematic approach to building com- 
plete ./V-representability conditions on the two-electron 
reduced density matrix (2-RDM) [28]. While an example 
of the new conditions was given previously, in the present 
paper we present further 7V-representability conditions 
on the 2-RDM that follow from the constructive solution. 
The conditions are in the form of a set of model Hamil- 
tonians with pairwise interactions whose trace against 
the 2-RDM must be nonnegative. The resulting condi- 
tions can be classified into an increasing hierarchy of con- 
straints, known as the (2, g)-positivity conditions where 
the first number p in the name indicates the highest p- 
RDM required to evaluate the condition (the 2-RDM in 
our case) and the second number q indicates the highest 
q-particle reduced density operators (g-RDOs) canceled 
by nonnegative linear combinations in the derivation of 
the condition. The (p, p)-positivity conditions are equiv- 
alent to the p-positivity conditions introduced earlier in 
Refs. [l(| EH Ell- We will use the two conventions in 
nomenclature interchangeably. 

In addition to the previously known Tl and T2 con- 
ditions [1, EE [I?], EM [HI, we derive a new class of 

(2.3) -positivity conditions. We also derive 3 classes of 

(2.4) -positivity conditions, 6 classes of (2,5)-positivity 
conditions, and 24 classes of (2,6)-positivity conditions. 
The conditions obtained can be divided into two general 
types: (i) lifting conditions, that is conditions which arise 
from lifting lower (2, g)-positivity conditions to higher 
(2, q + l)-positivity conditions and (ii) pure conditions, 
that is conditions which cannot be derived from a simple 
lifting of the lower conditions. All of the lifting conditions 
and the pure (2, g)-positivity conditions for q > 3 require 
that the expansion coefficients in the model Hamiltonians 
be tensor decomposed. Subsets of the A-representability 
conditions can be employed with previously known con- 



ditions for polynomially scaling calculations of ground- 
state energies and 2-RDMs of many-electron quantum 
systems in chemistry and physics. 

II. THEORY 

After the constructive solution of iV-representability is 
reviewed in section III Al it is employed in sections III CI 
and lll Dl to derive known and new A-representability con- 
ditions, respectively. The new constraints are organized 
into sections on (2,3)-, (2,4)-, (2,5)-, and (2,6)-positivity 
conditions. Two algorithms for implementing the condi- 
tions in a variational 2-RDM calculation are briefly dis- 
cussed in section Hi Bl 

A. Constructive solution 

The energy of an .ZV-electron quantum system in a sta- 
tionary state can be computed from the Hamiltonian 
traced against the state's density matrix 

E = Tt(H n D), (1) 

where the Hamiltonian operator is expressible in second 
quantization as 

H = Y. 2R M^ (2) 

ijkl 

in which the matrix 2 if is the reduced Hamiltonian oper- 
ator in a finite one-electron basis set [4(| and the indices 
label the members (orbitals) of the basis set. Because 
electrons are indistinguishable with pairwise interactions, 
the energy can also be universally written as a linear func- 
tional of only the 2-RDM 

E = Tt(H 2 D), (3) 

where the 2-RDM can be formally defined from integra- 
tion of the A-electron density matrix over all electrons 
save two 

2D= mN-Dj NDd3 dN (4) 

The expression of the energy as a functional of the 2- 
RDM suggests the tantalizing possibility of computing 
the ground-state energy of any electronic system as a 
functional of only the 2-RDM 0,11 EJ. Early calcu- 
lations by Coleman [4j, Tredgold [42[, and others, how- 
ever, showed that minimization of the energy as a 2-RDM 
functional produces unphysically low energies without 
additional constraints on the 2-RDM to ensure that it 
represents an TV-electron density matrix. In 1963 Cole- 
man called these constraints the A-representability con- 
ditions [4j. 

Building upon work by Garrod and Percus [Bj], Rum- 
mer in 1967 showed by the bipolar theorem j43| that there 
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exists a convex set (cone) of two-body operators { 2 Oi} 
whose trace against a potential 2-RDM will be nonnega- 
tive 

Tr( 2 O 2 L>)>0 (5) 

if and only if the 2-RDM is A-representable [(J. Hence, 
the set of two-body operators { 2 Oi} defines the set P^ 
of A-representable 2-RDMs. We say that the set { 2 Oi} 
is the polar of Pfj and denote it as P^ . Characterizing 
the set Pjj of 7V-representable 2-RDMs, therefore, would 
be complete if we could characterize its polar set P^ . 
Kummer's original result demonstrates the existence of 
the set Pff* , but it does not provide a prescription for 
constructing it. 

Recently, a constructive solution to the N- 
representability problem has been derived through 
the complete characterization of the polar set P^* |28| . 
In Ref. [28[ it is proven that the second-quantized 
representation of the operators { 2 Oi} in Pjy can be 
explicitly constructed as follows 

2 6 = wAdJ (6) 

i 

where Ci are polynomials in the creation and/or anni- 
hilation operators of degree less than or equal to r (the 
rank of the one-electron basis set) and u>i are nonnegative 
integer weights. The proof relies on the fact that P^ is 
contained within the set Pjy* of operators of degree < 2r 
whose trace against an A-electron density matrix must 
be nonnegative. Because the extreme elements (rays) of 
the convex cone P^* are readily expressed as [4J] 

QCj, (7) 

the extreme elements (rays) of P^* can be constructed 
from the conic combinations (or nonnegative linear com- 
binations) given in Eq. ©. The conic combinations, if 
divided by J2i w i i can be interpreted as convex combina- 
tions. Conic combinations are contained in P^* if and 
only if they cancel all three- and higher-body operators, 
that is polynomials in creation and annihilation opera- 
tors of degree greater than of equal to 6. 



in which the necessary set of operators (model Hamil- 
tonians) Oj, defining the boundary of the convex set of 
2-RDMs, must be determined itcrativcly. Given an ini- 
tial set of model-Hamiltonian constraints that bound the 
minimum energy, the three key steps of the algorithm are: 
(i) solving the linear program for the optimal 2-RDM, (ii) 
updating the set of model-Hamiltonian constraints in the 
linear program, and (iii) repeating steps (i) and (ii) un- 
til the 2-RDM is nonnegative in its trace with all model 
Hamiltonians explored in step (ii). In the second step, 
the trace of each model Hamiltonian with the 2-RDM is 
minimized by optimizing the Hamiltonian's parameters 
(expansion coefficients) , and if the final trace is negative, 
the model Hamiltonian with its optimized parameters is 
added to the constraints in Eq. (0). In practice, only a 
subset of model Hamiltonians from the constructive so- 
lution is employed. 

Some of the iV-representability constraints can be col- 
lected together as a single semidefinite constraint on the 
2-RDM. The generalization of a linear program to in- 
clude semidefinite constraints is known as a semidefinite 
program, and the solution of such a program is called 
semidefinite programming flH |46| . Efficient large-scale 
semidefinite programming algorithms have been devel- 
oped for the variational calculation of the 2-RDM 
HSJIMHIizifii- While the model Hamiltonians cor- 
responding to previously known A-representability con- 
ditions in section III CI can be expressed as semidefinite 
constraints, the model Hamiltonians corresponding to 
the new conditions in section Til Dl which use tensor de- 
compositions of the expansion coefficients in the C% oper- 
ators, cannot be written as traditional semidefinite con- 
straints. In practice, however, we can add these non- 
standard constraints to a semidefinite program contain- 
ing the standard semidefinite constraints by the three- 
step iterative procedure discussed above for the linear 
program. A main advantage of this second algorithm 
is that a large number of model Hamiltonians can be in- 
cluded by a single semidefinite constraint. A similar algo- 
rithm, to which we refer for further details, was proposed 
in Ref. [50j for imposing the T2 condition by recursively 
generated linear inequalities. 

C. Known conditions 



B. Practical implementation 

Before developing known and new 7V-representability 
conditions in sections III CI and III Dl respectively, in this 
section we briefly indicate their practical applications by 
sketching two algorithms for computing the ground-state 
2-RDM. Minimizing the ground-state energy as a func- 
tion of the 2-RDM constrained by these conditions can 
be formulated as a linear program 

minimize E = Tr(H 2 D) (8) 
such that Tr(dj 2 D) > for all j, (9) 



All previously known ./V-representability conditions are 
generated by the constructive solution. The most impor- 
tant representability conditions on the 2-RDM, derived 
by Coleman Q and Garrod and Percus [Bj], are the D, Q, 
and G conditions — also, known as the 2-positivity con- 
ditions [ll|. These conditions restrict the two-particle 
RDM 2 D, the two-hole RDM 2 Q, and the particle-hole 
RDM 2 G to be positive semidefinite, that is 

2 DhO (10) 
2 QhO (11) 
2 GhO, (12) 
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where the elements of the RDMs are given by 



u kl 



(^\a\a]aia k \^) 



2 q«h*Iw!4i*> 

2 G% = <*|a&-a}a fc |*) 



(13) 
(14) 
(15) 



and M y indicates that the matrix M is constrained 
to be positive semidefinite. Physically, these conditions 
correspond to constraining the probability distributions 
of two particles, two holes, as well as one particle and one 
hole to be nonnegative. The 2-positivity conditions are 
generated from the constructive solution by restricting 
the following three two-body operators from Eq. ^ to 
be nonnegative for all coefficients 6y 



2 6 D 


= 6 D & D 


(16) 


2 o Q 


= CqCq 


(17) 


2 d G 


= c G & G 


(18) 



C D 


= ^ 6 y oj a] 




ij 


Cq 


= ^2 did j 




ij 


Cg 


= ^bijdl&j 




ij 



where the Co, Cq, and Cg cover all polynomials in cre- 
ation and annihilation operators of degree two 

(19) 
(20) 
(21) 

Note that conic combinations are not present in these 
conditions because when the Ci operators are of degree 
2, the expectation values of the Oi operators only involve 
the 2-RDM 0. 

The other previously known ■/V-representability 
conditions— the Tl and T2 conditions [1, EH E3, El, Hlj — 
are part of the (2,3)-conditions that follow from the 
constructive solution. These semidefinite conditions 
on the 2-RDM are obtainable from conic combinations 
of three-particle metric matrices that cancel their 
dependence on the 3- RDM 0, 03] 



Tl : 
T2 



3 D + 3 QhO 
3 E + 3 FyQ. 



(22) 
(23) 



where in second quantization the matrix elements of 
these metric matrices are definable as 



,,„, - {^\d\dlald s d q d p \^) 



3 K% = (*|atata fe ata ? a p |*) 

3 Fpi k s = (*\a p a q ata k a]dl\V) 

3 Q'p J q k s = (*|a p a g a s a+atat|*) 



(24) 
(25) 
(26) 
(27) 



The four metric matrices 3 D, 3 E, 3 F, and 3 Q correspond 
to the probability distributions for three particles, two 
particles and a hole, one particle and two holes, and three 
holes, respectively [HE! 03]. Restricting the 3 D, 3 E, 



3 F, and 3 Q matrices to be positive semidefinite generates 
the 3-positivity conditions [Tl], 13 which depend on the 
3-RDM. While the Tl and T2 conditions are a subset of 
the 3-positivity conditions, they depend only upon the 
2-RDM because the 3-particle parts of 3 D and 3 Q (and 
3 E and 3 F) cancel upon addition [l7l[T^. For example, 
the matrix elements of Tl are given by 



rp-l ijk 



6 3jijk 
pqs 



18 1 D i p A 



• • (28) 



where P I is the p-particle identity matrix and A denotes 
the Grassmann wedge product [4(| HH . 

While the Tl condition is unique, three distinct forms 
of the T2 condition can be generated from rearranging 
the second-quantized operators in the definition of the 
F metric matrix relative to those in the 3 E metric ma- 
trix [1 81 ] . Consider the two variants of the 3 F matrix with 
the following matrix elements: 



3 j™* = <*|a p ata 3 ata fe aj|*) 
3 Fpl k s = (*|ata p a 3 atata fe |*). 



(29) 
(30) 



The 3-positivity condition 3 F y implies both 3 F y 
and 3 F y because reordering the creation and anni- 
hilation operators does not change the vector space cov- 
ered by the metric matrix. Changing the ordering of the 
second-quantized operators in the 3 F matrix relative to 
those in the 3 E matrix, however, does generate two ad- 
ditional T2 conditions 



T2 = 3 E 
T2 = 3 E 



f yo 
3 pyo. 



(31) 
(32) 



It was the T2 form of the T2 condition that was originally 
implemented by Zhao et al. [l5[ and Mazziotti [13 EH • 

The three T2 conditions are generated in the construc- 
tive solution by keeping the following two-body operators 
from Eq. ^ nonnegative 



(33) 
(34) 
(35) 



2 C>T2 


= c E c E - 


vc F c\ 


2 (1- 


= c E c E - 


V C F Cp 


2 o~ 

^T2 


= c E c\- 


^c F c F 



where 



Cf 



c F = 



Cp = 



ijk a i a j a k 
ijk 

^ b *ijk 
ijk 



a,ia d\ 



^ijk a i a k a 3 



Cz = 



5> - 

ijk 



ijk 



b*j k d\didj. 



(36) 
(37) 
(38) 
(39) 



The three T2 conditions can be combined into a single 
generalized T2 condition as shown in Refs. [HI, HH. The 
Tl condition is also produced in the constructive solution 
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by keeping the following two-body operator from Eq. 
nonnegative 



where 



2 d T1 = c D c D + c Q c Q 



ijk 

Cq = ^b* jk a i a ] a k . 

ijk 



(40) 

(41) 
(42) 



Because the second-quantized operators in Cd and Cq 
are anticommutative, there is only one Tl condition. Un- 
like the D, Q, and G conditions, both Tl and T2 con- 
ditions arise from the conic combination of a pair of 3- 
positive operators that cancels their dependence on the 
3- RDM. 



D. New conditions 

The constructive solution also produces new N- 
representability conditions on the 2-RDM (28|. In this 
section we will discuss the further conditions on the 2- 
RDM that emerge from conic combinations of three- 
, four-, five-, and six-particle operators in Eq. ([6]), 
which we denote as (2,3)-, (2,4)-, (2,5)-, and (2,6)- 
positivity conditions, respectively. All of the new N- 
representability conditions require a nonlinear factoriza- 
tion of the expansion coefficients to cancel the higher- 
particle operators. 



the indices on the tensor coefficients which is alphabeti- 
cal in both C(i, j, k) in Eq. (gBJ and C(k, in Eq. ([52]), 
Although the summation indices within C and its adjoint 
are distinct, we only show primes on the indices of the 
adjoint when the indices of the two operators appear in 
the same sum. Finally, for the ./V-representability condi- 
tions to be valid for real symmetric and general Hermi- 
tian RDMs, one-index tensors dk and dj, denote dk and 
d k , respectively. For multi-index tensors we employ the 
convention that the first subscript determines conjugacy, 
that is t>ij..rn — bij..m and o^j rn — b*j m . 

The first operator 0(i,j,k) arises from lifting the D 
condition through the insertion of a particle projection 
operator 



E 

fe,*;' 



d k d* k ,a\'a k > 



(47) 



while the second operator 0{i, j, k) arises from lifting the 
D condition through the insertion of a hole projection 
operator 



did 



k ak'0<ka k , 



(48) 



The nonnegativity of 0(i,j, k) and Q(i, j, k) gen erates a 
pair of lifting conditions discussed in Refs. [13j,|39(. While 
these two conditions depend not just on the 2-RDM but 
on parts of the 3-RDM, the sum of these two three-body 
operators produces a two-body operator 



2 o L1 = d(i,j,k) + 6(i,j,k). 



(49) 



1. (2,3)-positivity conditions 

In addition to the Tl and T2 conditions there exists 
a second class of (2,3)-positivity conditions that can be 
generated from lifting the 2-positivity conditions to the 
three-particle space and then canceling the three-particle 
operators. Consider the pair of three-body operators 



d(i,j,k) = c(i,j,k)6(i,j,k)i 
d(i,j,k) = c(i,j,k)c(i,j,h~y 



where 



C(i,j,k) 
C(i,j,k) 



ijk 

^hjd* k d\a\d k . 

ijk 



(43) 
(44) 

(45) 
(46) 



The notation for the operators 0(i, j, k) and C(i, j, k) in- 
cludes their internal summation indices to indicate suc- 
cinctly: (i) the ordering of the second-quantized opera- 
tors with indices i, j, and k, and (ii) the type of second- 
quantized operator with k denoting at, and k denoting dk- 
Note that the notation does not indicate the ordering of 



Because the two-body operator 2 Oli simplifies to the 
two-body operator 2 Ob in Eq. (|16[) . its nonnegativity 
regenerates the D condition. With a generalization of this 
lifting process, however, we can generate (2,3)-positivity 
conditions that are distinct from the known conditions. 

We can generalize the lifting process by inserting the 
creation operator and the annihilation operator respon- 
sible for lifting at non-adjacent positions. For example, 
consider the pair of three-body operators 



where 



6(k,ij) = c(k,i,j)c(k,i,tf 
d(k,i,j) = c(k,ij)6(k,ijy 

C(k, = K d ka\d\d] 

ijk 

C(k,i,j) = Yhjdld k d\d). 

ijk 



(50) 
(51) 

(52) 
(53) 



In 0(k,i,j) the creation operator d k in C and the an- 
nihilation operator dk' in the adjoint of C , which per- 
form the lifting of the D condition, are separated from 
each other by four second-quantized operators; similarly, 
in 0(k,i,j) the creation and annihilation operators, dk 
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TABLE I. The (2,3)-positivity conditions can be derived from conic (linear nonnegative) combinations of the (3,3)-positivity 
conditions that cancel the 3-particle operators. 



Class 


Type 


Representative Condition 


C Definition 


1 

2 


Lifted (2,2) 
Pure (2,3) 


Tr((0(k,i,j) + 0(k,i,j)) 2 D) >0 
Tr((6(i,j,k) + 6(i,j,k)) 2 D)>0 


Eq. (HD 
Eq. (HQ 



and a k , respectively, are separated from each other by 
four second-quantized operators. Because the compo- 
nents of the projectors are separated, the nonnegativ- 
ity of 0(k,i,j) and 0(k,i,j) generates a pair of gener- 
alized lifting conditions that extend those discussed in 
Refs. [II Eg. 

While individually 0(k,i,j) and 0(k,i,j) depend on 
three-particle operators, their sum generates a two-body 
operator 

2 d L2 = 6(k,ij) + d(k,i,j). (54) 

Unlike 2 Oli, the nonnegativity of the lifted operator 
2 Ol2 is not necessarily implied by the D, Q, G, Tl, and 
T2 conditions. Importantly, 2 Ol2 does not simply re- 
arrange to 2 Oli because the creation and annihilation 
operators are non-commutative. Based on the possible 
orderings of the fundamental second-quantized opera- 
tors, there are nine distinct ways to lift the D condition 
while canceling the resulting three-particle operators and 
hence, nine distinct lifting conditions from the D condi- 
tion. Similarly, there are nine distinct (2,3)-positivity 
conditions from lifting the Q condition and nine from 
lifting the G condition. Three of these 27 lifting con- 
ditions reduce to the D, Q, and G conditions, respec- 
tively while the other conditions are distinct because the 
second-quantized operators in quantum mechanics form 
a non-commutative algebra. 

Table I summarizes the (2,3)-positivity conditions by 
giving a representative condition from each of the two 
classes: (i) the lifting conditions and (ii) the pure con- 
ditions. While the lifting conditions arise from lifting 
the 2-positivity conditions, the pure conditions cannot 
be obtained from lifting any of the lower conditions. Ta- 
ble I gives nonnegativity of 2 Ol2 and the Tl condition 
as representative conditions of the lifting and pure (2,3)- 
positivity conditions, respectively. All of the other (2,3)- 
conditions can be obtained from these representative con- 
ditions through two processes, (i) switching of the second- 
quantized operators in the C(i, j, k) between creators and 
annihilators and (ii) reordering of the second-quantized 
operators in the C(i,j, k). 

Switching the second-quantized operators with index j 
in the L2 condition of Eq. (|54l) . for example, generates a 
lifted G condition 

2 6 L3 = d(k,ij) + 6(k,i,j). (55) 

Note that switching the second-quantized operators with 
index k in Eq. (|54p simply regenerates the same condi- 



tion while switching the second-quantized operators as- 
sociated with indices i and j generates lifted G and Q 
conditions from the lifted D condition. Reordering of 
the second-quantized operators in the L2 condition of 
Eq. (|54l) . by contrast, produces the other 9 lifted D condi- 
tions; for example, reordering L2 yields the LI condition 
in Eq. (j49|) . Similarly, for the pure (2,3)-positivity con- 
ditions switching of the second-quantized operators with 
index k in Tl produces the T2 condition in Eq. (l23l) . The 
other two distinct T2 conditions, T2 and T2, in Eqs. (|31~j) 
and (1321) are generated not by switching but by reorder- 
ing the second-quantized operators in the T2 condition 
of Eq. 423). 

2. (2,4)-posittvity conditions 

The (2,4)-positivity conditions, arising from consider- 
ing all Ci operators of degree less than or equal to four 
in Eq. ((6]), consist of two classes of lifting conditions and 
one class of pure conditions, which are summarized in 
Table II. The two classes of lifting conditions are gener- 
ated from lifting the two classes of (2,3)-positivity condi- 
tions. As in the previous section, the generalized lifting 
is performed by (i) inserting a creation operator into each 
Ci operator contributing to the condition, (ii) converting 
the inserted creation operator into an annihilation op- 
erator in the operator produced from step (i), and (iii) 
adding the two lifted operators from steps (i) and (ii) 
together to produce a two-particle operator. The non- 
negativity of the resulting two-particle operator gener- 
ates a lifting (2,4)-positivity condition. Representative 
lifting conditions for both classes are shown in Table II. 
The C operators in the first and second classes of lifted 
(2,3)-positivity conditions and the pure (2,4)-positivity 
condition are given by 

(7(7, k, i, j) = ^2 bijd k eia\a\a\a] (56) 

ijkl 

C(l, i, j, k) = ^2 hjkdia\a\a)a\ (57) 

ijkl 

C(i,j,k,l) = '^b i d ] e k fiaJ i a) j d\a] (58) 

ijkl 

where d;, e^, fc, bij, b ijkj and aj become b*, d*, e*, /*, 
b*j, b*j k , and di when i = i. The rank of the largest ten- 
sor changes from three in Eq. (|5"6")) to one in Eq. (1551) to 
effect the cancelation of the 3- and 4-RDOs in the com- 
binations of operators in Table II. Fusing the tensors bi 
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TABLE II. The (2,4)-positivity conditions can be derived from conic (linear nonnegative) combinations of the (4,4)-positivity 
conditions that cancel the 3- and 4-particle operators. 



Class Type 


Representative Condition 


C Definition 


1 Lifted (2,2) 

2 Lifted (2,3) 

3 Pure (2,4) Tr((30(i 


Tr((6(I, k, i,j) + 6(1, k, + 6(1, k, + 6(1, k, i, j)) 2 D) > 
Tr((6(Z, i, j, k) + 6(l,i,j,k) + 6(1, i,j,k) + 6(1, i, j, k)) 2 D) > 

, j, k, i) + 6(i, j, k, i) + 6(i,j, k, i) + 6(iJ, k, i) + 6(1 j, k, i) + 6(1 j, k, I)) 2 D) > o 


Eq. ([56} 
Eq. J57J 
Eq. (JSHJ 



and dj into a single rank- two tensor by in Eq. ff58|) . for 
example, would cause the operator combinations in Ta- 
ble II to depend on the 3- and 4-RDOs. Additional (2,4)- 
positivity conditions can be generated from the represen- 
tative conditions through a combination of switching and 
reordering of the creation and annihilation operators. 

The pure (2,4)-positivity conditions, presented in 
Ref. |2g], depend upon only the 2-RDM through conic 
combinations that cancel the 3- and 4-RDMs. As in the 
(2,3)-conditions, the cancelations depend upon the conic 
combination of pairs of operators that differ from each 
other by an odd number of switchings — exchanges of cre- 
ators and annihilators. Generating an extreme condition 
on the 2-RDM requires that we consider the minimum 
number of conic combinations that effect the cancelation 
of the higher RDMs. Each pure (2,4)-positivity condi- 
tion involves the conic combination of eight four-particle 
operators by Eq. ©. These eight four-particle operators 
can be grouped into the four pairs that depend upon only 
three-particle operators: 

d(i,j,k,i) + d(i,],k,i) (59) 
d(i,],k,i) + d(i,j,k,i) (60) 
d(i,j,k,i) + 6(i,j,k,i) (6i) 

6(i,j,k,l)+6(i,j,k,l) (62) 

The operators in the first pair differ from each other by 
the switching of three creation and annihilation opera- 
tors while the operators in the other three pairs differ 
from each other by the switching of one creation operator 
and one annihilation operator. Rearranging the second- 
quantized operators in the four pairings into normal or- 
der with creators to the left of the annihilators generates 
expressions involving the sum of 9, 5, 3, and 1 3-RDOs, 
respectively. Upon summation the 9 3-RDOs from the 
one pairing with three switchings cancels with the 5, 3, 
and 1 3-RDOs from the three pairings with one switch- 
ing, and hence the final operator depends upon only the 
2-RDO. 

Other pure (2,4)-positivity conditions can be gener- 
ated from the representative condition through switch- 
ing and reordering of the second-quantized operators. 
To maintain the cancelation of the 3- and 4-RDOs, we 
must perform the same switching of creation and anni- 
hilation operators in each operator C(i,j, k, I) contribut- 
ing to the condition. Because each fundamental second- 
quantized operator can be either a creation or an an- 
nihilation operator, there are 2 4 or 16 conditions from 



switching. Eight of these conditions can be generated 
from the other 8 conditions by switching all creation and 
annihilation operators by particle-hole symmetry. In the 
limit that the expansion coefficients bi , dj , e k , and /; be- 
come orthogonal unit vectors, these 16 conditions reduce 
to the 16 conditions in (2,4)-class of the classical (or diag- 
onal) A^-representability problem 0, [3(1 H3| • The quan- 
tum mechanical formulation of these conditions, however, 
is much more general because the expansion coefficients 
need not be orthogonal. When the expansion coefficients 
are non-orthogonal, the creation and annihilation opera- 
tors become non-commutative operators, and hence, the 
conditions depend upon their ordering. 

In the quantum case additional conditions can be gen- 
erated from each of the 16 conditions by reordering the 
creation and annihilation operators while preserving the 
cancelation of the 3- and 4-RDOs. These additional con- 
ditions are related to the original 16 conditions as the 
generalized T2 conditions are related to the T2 condi- 
tion in the (2,3)-positivity conditions. Table III presents 
the representative pure (2,4)-positivity condition as well 
as three other conditions generated from its reordering. 
Each of these four conditions differs from the others by 
a few terms involving the 2-RDM. For example, the first 
and second conditions differ by only one term 

92 = 5i + I a(3 £ he k 2 D% b*,/;, > 0. (63) 

\ ikyi' V J 

where 

a = £>e* (64) 

3 

P = Y,f3 d *3 ( 65 ) 

3 

and 5? selects the real part of the expression. When this 
term is negative, inequality g2 is stronger than g\, but 
when this term is positive, inequality g\ is stronger than 
92- In the classical case, where the expansion coefficients 
are orthogonal, these two conditions are equivalent be- 
cause both a and /3 are zero, and hence, this additional 
term vanishes. 

3. (2,5)-positivity conditions 

The (2,5)-positivity conditions are generated from con- 
sidering all Ci operators of degree less than or equal to 
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TABLE III. The representative pure (2,4)-positivity condition g\ > as well as three other conditions generated from its 
reordering, Q2 > 0, gs > 0, and gi > 0, are presented. Unlike the situation in the classical limit, in the quantum case additional 
conditions can be generated from each of the 16 conditions obtained from switching by reordering the creation and annihilation 
operators while preserving the cancelation of the 3- and 4-particle operators. 



Condition C Definition 

9l ( 2 D) = Tr((30(i, j, k, T) + 0{i, j, k, T) + 6(i, j, k, I) + 0(i, j, k, I) + 0(i, j, kj) + 0(i,j,k,l)) 2 D) > Eq. J58J 

g 2 { 2 D) = Tr((30(i, j, k, I) + 6(i,j, k, I) + 6(i, k,j, I) + 6(k,j, i, I) + 6(j, i, k, I) + 0(i,j, k, lj) 2 D) > Eq. (GU 

g 3 ( 2 D) = Tr((36(i, j, k, I) + 6(i, j, l,k) + 6(i,l,k, j) + 6(k,j, i, I) + 6(j, i, k, I) + 6(i, j,k,l)) 2 D)>0 Eq. J58J 

g 4 ( 2 D) = Tr((36(i, j, k, I) + 6(i, j, k,l) + 6(i, k, j, I) + 6(k, j, i, I) + 6(j, i, k,l) + 6{i,j,k,l)) 2 D) > Eq. J58} 



TABLE IV. The (2,5)-positivity conditions can be derived from conic (linear nonnegative) combinations of the (5,5)-positivity 
conditions that cancel the 3, 4-, and 5-particle operators. 



Class 



Type 



Representative Condition 



C Definition 



1 Lifted (2,2) 

2 Lifted (2,3) 

3 Lifted (2,4) 

4 Pure (2,5) 

5 Pure (2,5) 



6 Pure (2,5) 



Tr((0(m,l,k,i,j) + 0{m,l, k,i, j) + 0(fh,l,k,i,j) + 0(rh,l,k,i,j) Eq. JMJ 

+ 6{m,l,k,i,j) + 6(m,l,k,i,j) + 6(fh,l,k,i,j) + 0{m,l,k, i, j)) 2 D) > 
Tr((6(m,Z, i,j,k) + 6{m,l,i,],k) + 6{m, 1,1, ], k) + 0{rh,l,i,j,k) Eq. dHZJ 

+ d(m,l,i,j,k) + 6{m,l,i,j,k) + 6(rh,l,i,j,k) + 6(rh,l,i, j, k)) 2 D) > 

Tr((30{m,i,j,k,l) + 6(m,i,j,k,l) + 6(m,i,j, k,l) Eq. J68| 

+ 0(m, i, j, k, I) + 0(m, i, j, k, I) + 0(m, i,j,k,l) 
+ 3(5(m, i,j,k,l) + 0(fh,i,j,k,l) + 0(fh,i,j,k,l) 
+ 6(m,i,j,k,l) + 6(m,i,j,k,l) + 6(rh,i,j ,k,l)) 2 D) > 
Tr((3d{i,j,k,l,m) + d(i,j,k,l,rh) + d(i,j,k,T,m) Eq. (|69) 

+ 0(i,j, k, I, m) + (D(i,j, k, I, m) + 0(i,j, k, I, m) 
+ 30(i,j, k, I, rh) + 0(i,j, k, I, m) + 0(i,j, k, I, m) 
+ 6(i,j,k,l,m) + d(t,j,k,l,fh) + d(i,],k,l,rh)) 2 D) > 
Tr((60{i,j,k,l,m) + 3d\i,j,k,l,fh) + 30{i,j,k,l,m) Eq. (|69) 

+ 30(i,j, k, I, m) + 30(i,j, k, I, m) + 30(i,j, k, I, m) 
+ 0(i,j, k, I, m) + 0(i,j, k, I, m) + 0(i,j, k, I, m) 
+ 0(i,j, k,l,fn) + 0(i,j, k, I, m) + 0(i,j, k,l,m) 
+ 0(i, j, k,l,m) + 0(i,j, k, I, m) + <D(i,j, k, I, m) 
+d{i,],k,l,m) + 6{i,j,k,l,rh)) 2 D) > 
Tr((60(i,j,k,l,m) + 30(i,j,k,l,rh) + 30(i,j,k,l,m) Eq. (|69} 

+ 30(i,j, k,l,m) + 30(i,j, k, I, m) + 0(i,j, k,l,m) 
+ 0(i,j, k, I, m) + 0(i,j, k, I, fh) + 0(i,j, k, I, m) 
+ 0(i,j,k,l, fh) + 0(i,j,k,l, m) + 0(i, j ,k,l, m) 
+ 0(i, j, k,l,fh) + 0(i,j, k, I, m) + 0(i,j, k, I, m) 
+6{i,j, k, I, m) + 30(i,j, k, I, fh)) 2 D) > 



five in Eq. ((6]). These conditions consist of three classes 
of lifting conditions and three classes of pure conditions, 
which are given in Table IV. The lifting conditions arise 
from lifting the three different classes of (2,4)-positivity 
conditions. The C operators of the first, second, and 
third classes of lifting conditions are given by 



C(m,l,k,i,j) = ^ hjdkeifmat n a}al.ala' j (66) 

ijklm 

C(m,l,i,j,k)= ^ kjkdiemalnaja^atal (67) 

ijklm 

C(m,i,j,k,l) = ^ b l d j e k fig m ai l a\a'lala}, (68) 

ijklm 



respectively, and the C operators of the three classes of 
pure conditions are given by 

C(i,j,k,l,m)= 2J bidjekfigrnolalalajaln (69) 

ijklm 

where the b i} di, e,, fi, g l , bij, b ijk , and a\ become b*, 
d*, e*, /*, g*, b*j, b* jk , and a, when i = i. Note that the 
operators in Eqs. (ffjg|) and (|69|) are not equivalent after 
switching. Switching of creators to annihilators in the 
C operators in the representative conditions in Table IV 
produces 16, 32, and 32 conditions in the pure classes 
4, 5, and 6, respectively. Class 4 has fewer conditions 
because its conditions, unlike those in classes 5 and 6, 
possess particle-hole symmetry. Particle-hole symmetry 
is present in all of the pure (2,3)-positivity conditions 
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and none of the pure (2,4)-conditions. Additional condi- 
tions can be generated from the representative conditions 
through reordering of the creation and annihilation oper- 
ators. Like the (2,3)- and (2,4)-positivity conditions, the 
(2,5)-conditions generate all of the classical (diagonal) 
7V-representability conditions when the expansion coeffi- 
cients bi, dj, efe, fi, and g m are chosen to be orthogonal 
unit vectors. 



4- (2,6)-positivity conditions 

As with the (2,g)-positivity conditions for q < 5, the 
(2,6)-positivity conditions are generated from Eq. © by 
considering all operators of degree less than or equal to 
six. Six classes of lifting (2,6)-positivity conditions arise 
from lifting the six classes of (2,5)-positivity conditions. 
While not shown explicitly, the representative conditions 
can be readily constructed from the conditions in Ta- 
ble IV. There are also 18 classes of pure (2,6)-positivity 
conditions. The C operators of these 18 conditions are 
given by 

C{ijklmn) = hdjekfigmhnald^dlaja^a^ (70) 

ijklmn 

where the bi, di, &i, fi, Qi, hi, and a\ become b*, d*, 
e*, /*, g*, h*, and dj when i = i. Table V provides a 
representative operator for each of the 18 classes. Each 
representative operator arises from the conic combination 
of potentially 2 6 (or 64) six-particle operators, which are 
distinguished from each other by the switching between 
creation and annihilation operators. These 64 operators 
are grouped in 32 particle-hole pairs given in the rows 
of Table V. For each of the 18 representative conditions, 
the nonnegative integer weights a and /3 of the operators 
in each pair are reported. The conic combination of all 
32 pairs with the weights in the x th column generates a 
representative operator for class x. The operator for each 
class depends only on the 2-RDO with the dependence 
on the 3-, 4-, 5-, and 6-RDOs canceling through the conic 
combination. The trace of each representative operator 
against the 2-RDM generates a representative condition 
on the 2-RDM. Additional (2,6)-positivity conditions can 
be generated from the representative conditions through 
a combination of switching and reordering of the creation 
and annihilation operators. From the particle-hole pair- 
ing it is easy to observe that only one class of the (2,6)- 
conditions — class 4 — has particle-hole symmetry, that is 
a = (3 in all pairs. 

The (2,6)-positivity conditions yield all classes of the 
classical (diagonal) A^-representability conditions when 
the expansion coefficients bi, dj, efe, fi, g m , and h n are 
chosen to be orthogonal unit vectors. Classically, all 
classes of (2,q)-conditions for o < 5 are in the form of 
hypermetric inequalities [36l [37| . When q = 6, how- 
ever, new classes of classical Af-representability condi- 
tions emerge [1, [3(1 H3, HH • In the classical limit, the first 



6 classes of pure (2,6)-positivity conditions in Table V re- 
duce to hypermetric inequalities while the remaining 12 
can be grouped into cycle, parachute, and Grishukhin 
inequalities [52j. 



III. DISCUSSION AND CONCLUSIONS 

Both new and known A^-representability conditions on 
the 2-RDM have been derived from the constructive so- 
lution to the iV-representability problem 28]. In addi- 
tion to all of the previously known conditions, we gener- 
ate new (2,3)-, (2,4)-, (2,5), and (2,6)-conditions where 
the first number p in each pair indicates the highest p- 
RDM required to evaluate the condition (the 2-RDM in 
our case) and the second number q indicates the highest 
RDMs canceled by conic (linear nonnegative) combina- 
tions in the derivation of the condition. There are two 
classes of (2,3)-conditions: (i) lifting conditions that are 
derivable from lifting the D, Q, and G (2-positivity) con- 
ditions to the three-particle space, and (ii) pure condi- 
tions that are not derivable from lifting and hence, are 
without precedent in the 2-positivity conditions. The 
(2,4)-conditions have two classes of lifting conditions and 
one class of pure conditions, the (2,5)-conditions have 
three classes of lifting conditions and three classes of pure 
conditions, and the (2,6)-conditions have six classes of 
lifting conditions and eighteen classes of pure conditions. 
A similar procedure of using conic combinations to can- 
cel operators higher than two-body can be followed for 
deriving the (2, g)-conditions for q > 6. 

The classical (diagonal) A^-representability condi- 
tions 0, are constraints on the two-electron re- 
duced density function, the diagonal part of the 2-RDM, 
to ensure that it represents an iV-electron density func- 
tion. A solution to the diagonal problem was devel- 
oped in the context of both the Boole 0-1 programming 
problem and the maximum cut problem of graph the- 
ory [37l [53j . The recent constructive solution of the N- 
representability problem for fermionic density matrices 
extends the classical solution to the more general quan- 
tum case. All of the quantum conditions can be cast 
in the form of restricting the trace of two-body oper- 
ators (model Hamiltonians) against the 2-RDM to be 
nonnegative. In the limit that all tensors in the model 
Hamiltonians are decomposed into products of orthog- 
onal rank-one (one-index) tensors, the quantum condi- 
tions reduce to the classical (diagonal) conditions for 
all unitary transformations of the one-electron basis set. 
The quantum (2,6)-conditions presented here reduce in 
the classical limit to the complete set of classical (2,6)- 
conditions 13(1 which were shown to be complete by 
Grishukhin [52j. 

A significant difference between the classical and quan- 
tum conditions is the orthogonality (classical) or non- 
orthogonality (quantum) of the rank-one tensors. Con- 
sequently, in the classical case the creation and anni- 
hilation operators form a commutative algebra while in 
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TABLE V. The (2,6)-positivity conditions can be derived from conic combinations of the (6,6)-positivity conditions that cancel 
the 3, 4-, 5-, and 6-particle operators. There are 6 classes of lifting conditions (not shown) and 18 classes of pure conditions 
(shown). The table provides a representative operator for each of the 18 classes. The conic combination of all 32 pairs with 
the weights in the x th column generates a representative operator for class x. The trace of each representative operator against 
the 2-RDM generates a representative condition on the 2-RDM. 



Weights (a/ /3) 



Operators 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 18 


aO(ijklmn] 


+ /30(ijklfhn) 


20/2 


12/6 20/6 12/12 20/12 20/12 9/6 


5/3 


9/5 


9/3 


5/2 


14/9 14/6 14/3 


6/3 


6/5 


1/6 3/3 


aO(ijklmn] 


+ f30(ijklrhn) 


12/0 


6/2 


12/2 


6/6 


12/6 


12/6 


5/3 


2/1 


3/1 


5/1 


3/1 


9/5 


9/3 


9/1 


5/2 


5/3 


0/3 2/1 


aO(ijklrhn] 


+ f30(ijklmn) 


12/0 


6/2 


12/2 


6/6 


12/6 


12/6 


6/3 


3/1 


6/3 


6/1 


3/1 


9/5 


9/3 


9/1 


3/1 


3/3 


0/3 1/1 


aO(ijklrhn] 


+ (30(ijklmn) 


6/0 


2/0 


6/0 


2/2 


6/2 


6/2 


3/1 


1/0 


1/0 


3/0 


1/0 


5/2 


5/1 


5/0 


3/1 


2/1 


0/1 1/0 


aO{ijklmn] 


+ /30(ijklfhn) 


12/0 


6/2 


12/2 


6/6 


12/6 


12/6 


5/3 


2/1 


5/2 


5/1 


3/1 


9/5 


9/3 


9/1 


3/1 


3/2 


0/3 3/2 


aO(ijklmn] 


+ /30(ijklrhn) 


6/0 


2/0 


6/0 


2/2 


6/2 


6/2 


2/1 


0/0 


1/0 


2/0 


1/0 


5/2 


5/1 


5/0 


2/0 


3/1 


0/1 3/1 


aO(ijklmn] 


+ /30(ijklmn) 


6/0 


2/0 


6/0 


2/2 


6/2 


6/2 


3/1 


1/0 


3/1 


3/0 


2/1 


5/2 


5/1 


5/0 


1/0 


1/1 


0/1 2/1 


aO(ijklmn] 


+ fiO(ijklmn) 


2/2 


0/0 


2/0 


0/0 


2/0 


2/0 


1/0 


0/0 


0/0 


1/0 


0/0 


2/0 


2/0 


2/0 


1/0 


1/0 


1/0 3/1 


aO(ijklmn] 


+ /30(ijklrhn) 


12/0 


6/2 


12/2 


6/6 


12/6 


12/6 


3/2 


3/2 


6/3 


6/2 


3/1 


6/3 


9/3 


6/0 


3/1 


3/2 


1/5 2/3 


aO(ijklmn] 


+ /3(J(ijklrhn) 


6/0 


2/0 


6/0 


2/2 


6/2 


6/2 


1/1 


1/1 


2/1 


3/1 


2/1 


3/1 


5/1 


3/0 


3/1 


3/1 


0/2 1/1 


aO(ijklrhn] 


+ /3(J(ijklmn) 


6/0 


2/0 


6/0 


2/2 


6/2 


6/2 


1/0 


1/0 


3/1 


3/0 


1/0 


3/1 


5/1 


3/0 


1/0 


1/1 


0/2 0/1 


aO(ijklrhn] 


+ f30(ijklmn) 


2/2 


0/0 


2/0 


0/0 


2/0 


2/0 


0/0 


0/0 


0/0 


1/0 


0/0 


1/0 


2/0 


1/1 


2/1 


1/0 


0/0 0/0 


aO(ijklmn] 


+ f30(ijklmn) 


6/0 


2/0 


6/0 


2/2 


6/2 


6/2 


1/1 


1/1 


3/1 


3/1 


2/1 


3/1 


5/1 


3/0 


1/0 


1/0 


1/3 1/1 


aO(ijklmn] 


+ /30(ijklfhn) 


2/2 


0/0 


2/0 


0/0 


2/0 


2/0 


0/1 


0/1 


1/1 


1/1 


1/1 


1/0 


2/0 


1/1 


1/0 


2/0 


i/i i/o 


aO(ijklfhn] 


+ f30(ijklmn) 


2/2 


0/0 


2/0 


0/0 


2/0 


2/0 


0/0 


0/0 


1/0 


1/0 


1/1 


1/0 


2/0 


1/1 


0/0 


0/0 


1/1 0/0 


aO(ijklmn] 


+ /30(ijklrhn) 


12/0 


6/2 


12/2 


6/6 


6/2 


12/6 


3/1 


3/1 


5/2 


5/1 


2/0 


3/1 


3/0 


9/1 


1/0 


1/1 


1/5 1/1 


aO(ijklmn] 


+ /30(ijklfhn) 


6/0 


2/0 


6/0 


2/2 


2/0 


6/2 


1/0 


1/0 


1/0 
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the quantum case they form a non-commutative algebra. 
The non-orthogonality leads to active TV-representability 
conditions on the 2-RDM that lack a classical analogue. 
For example, all classes of lifting conditions that we pre- 
sented are inactive in the classical limit. Because the 
creation and annihilation operators commute, each class 
of classical (2, g)-lifting conditions reduces to a class of 
classical (2, p) -pure conditions where p < q. Further- 
more, typically more than one pure quantum condition 
reduces to each classical condition in the classical limit. 
Table III shows four pure (2,4)-conditions that reduce to 
the same classical condition. These quantum conditions 
differ only in the ordering of the creation and annihilation 
operators — a difference that disappears in the classical, 
commutative limit. 

The conic combination of the extreme two-body oper- 



ators in the iV-representability conditions forms a con- 
vex set (cone) of model Hamiltonians for which the N- 
representability conditions are exact. From the perspec- 
tive of quantum information the computational complex- 
ity of enforcing all iV-representability conditions on the 2- 
RDM can be shown to be non-deterministic polynomial- 
time (NP) complete, meaning that in the worst-case 
scenario enforcing exact 7V-representability scales non- 
polynomially with system size. Despite this complexity, 
however, many realistic quantum systems are much more 
tractable than the worst-case scenario implies. For exam- 
ple, the 2-positivity conditions, particularly the G con- 
dition, are exact for pairing Hamiltonians whose ground 
states are antisymmetrized geminal power wavefunctions. 
Such pairing Hamiltonians have been employed to model 
the Cooper pairing and long-range order associated with 
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superconductivity. For any strength of interaction the 
ground-state energy for this class of Hamiltonians can be 
computed in polynomial time. 

More generally, for fixed q the (2, (^-positivity condi- 
tions, which contain the lower positivity conditions, cover 
a large class of model Hamiltonians whose ground states 
are computable in polynomial time — in a time that scales 
polynomially with system size. Even when the Hamilto- 
nian of interest is not rigorously contained in this class, 
the associated ./V-representability conditions, which in- 
trinsically are not constrained by the approximations 
of perturbation theory, may produce an accurate lower 
bound on the ground-state energy. Computational expe- 
rience with the variational calculation of the 2-RDM in 
atoms and molecules 0, Q M, H M, US E3 shows 
that sufficiently accurate lower-bound ground-state ener- 
gies are often produced with (2, g)-positivity conditions 
where q < 3. 

The practical implementation of the variational 2- 
RDM method requires that the energy be minimized 
as a functional of the 2-RDM constrained by its N- 
representability conditions. Both the 2-positivity con- 
ditions and the Tl and T2 conditions can be expressed 
as positive semidefinite constraints (also known as linear 
matrix inequalities) in which metric matrices are con- 
strained to be positive semidefinite. These constraints 
on the 2-RDM can be imposed during the minimization 
of the ground-state energy through a genre of constrained 
optimization known as semidefinite programming 
Es HI HI H US-SI- The remaining (2, g)-positivity 
conditions, however, cannot be expressed as a traditional 
semidefinite constraint because the coefficients in the Ci 
operators must be tensor decomposed to remove the de- 
pendence of the constraints on the higher RDMs. Prac- 
tically, as described in section HlBl these constraints can 
be added to the semidefinite program through recursively 
generated linear inequalities, similar to those described 
in Ref. |o| for T2. 

The constructive solution of iV-representability estab- 
lishes 2-RDM theory as a fundamental theory for many- 
particle quantum mechanics for particles with pairwisc 
interactions. Lower bounds on the ground-state energy 
can be computed and improved systematically within 
the theory. While not all of the 2-RDM conditions will 



be imposed in practical calculations, a complete knowl- 
edge of the conditions — their form and function — can 
be invaluable in devising and testing approximate N- 
representability conditions for different types of quantum 
systems and interactions. Like Feynman diagrams the 
positivity conditions represent different physical interac- 
tions of the electrons. Adding positivity conditions to the 
2-RDM calculation expands the class of exactly describ- 
able model Hamiltonians. Just as classes of Feynman dia- 
grams differ in importance according to the nature of the 
interaction, for a given system some positivity conditions 
will be significantly more important than others. For ex- 
ample, both the G and T2 conditions have proven to 
be especially important in calculations of many-electron 
atoms and molecules [3, [TBI, [ljj while the Tl condition 
has rarely been of any significance. Similar evaluations 
must be performed in a variety of many-electron quantum 
systems for the conditions resulting from the constructive 
solution. 

Previous variational 2-RDM computations on metal- 
lic hydrogen chains [5o| , polyaromatic hydrocarbons [25L 
[2"o| . and firefly luciferin [54| show that they can capture 
strong, mult i- reference correlation effects for which ap- 
propriate ansatze for the wavefunction are difficult to 
construct. With a suitable choice of 7V-representability 
conditions, therefore, strong electron correlation effects 
can be computed at a computational cost that scales 
polynomially with the system size. Although the ex- 
ploration of the conditions following from the construc- 
tive solution is still in its earliest stages, a 2-RDM- 
based theory with systematically improvable accuracy 
promises fresh theoretical and computational possibilities 
for treating strong correlation in quantum many-electron 
systems. 
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